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Abstract: 

We propose event cuts in deep inelastic scattering, suitable for an examination of quark 
fragmentation universality. We compare the current hemisphere of the Breit frame with 
hemispheres in e + e~ annihilation events, using Monte Carlo simulations. The agreement 
between the two processes is improved after the suggested events cuts. A method to study 
the scale evolution in quark hemispheres using data from fixed energy e + e~ experiments is 
presented. This makes it possible to use the high statistics from LEP1 also at scales below 
the Z° mass. We also discuss observables which are sensitive to the dynamics of regions 
closer to the remnant. The observables probes the relatively clean region on the current 
side of the hardest emission in the event, and can be used to distinguish between different 
assumptions about remnant effects and mechanisms for the parton evolution. 
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1 Introduction 



In deep inelastic scattering (DIS), the current hemisphere of the Breit frame is expected to 
be very similar to one hemisphere in an e + e~ experiment. This expectation relies on the 
fundamental assumption of quark fragmentation universality. However, a set of features in 
DIS introduce corrections to this assumption. 

QCD radiation can give rise to high-p^ emissions, where Q or higher. These emis- 
sions have no correspondance in an e + e~ event, where the kinematical constraint for gluon 
emission is p± < a/s/2. A high-p^ emission in DIS significantly changes the kinematical 
properties of the struck quark, and a high-pj_ event sometimes manifests itself with a com- 
pletely empty current region ||T|. If the rate of such events is getting non-negligeble, they 
may introduce an uncertainty to the interpretation of the data. Furthermore, problems arise 
even if the current region is not empty, but merely depopulated. It is therefore of interest 
to find a way to exclude high-pj_ events using another signal than an empty current Breit 
hemisphere. One candidate method is to reconstruct the highest p_i_-scale of the event using 
jets, and we will in this paper discuss that approach. 

The flavour composition in e + e~ and DIS is not exactly the same. When excluding high- 
p± events from the DIS analysis, the boson-gluon fusion channel for heavy quark production 
is supressed. This implies a lower heavy quark rate in the studied DIS sample, as compared 
to e + e~ data at corresponding energies. A uds enriched e + e~ data sample with high statistics 
is available from the Z° peak. We will in this paper discuss a method to study properties of 
e + e~ quark hemispheres at different scales, using data from the fixed energy Z° experiments. 

At very low squared momentum transfer, when Q 2 is of the same order as the squared 
p-meson mass m 2 , the vector meson component of the photon affects the current region 
properties. We will restrict the discussion to Q 2 > 4GeV 2 ~ 8m 2 . In this kinematical range 
the vector meson component of the probe is suppressed, and it will not be discussed further 
in this paper. 

Comparisons of DIS and e + e~ data support quark fragmentation universality ||. Some 
differences are however observed, especially at low Q 2 . These can in general be qualitatively 
explained from the features discussed above, but data on e.g. strangeness rates in DIS are 
less well understood 0. 

In this paper, we investigate how the similarity between e + e~ and current Breit hemi- 
spheres can be improved with proper event cuts. The effects of event cuts are examined 
with MC simulations. We also discuss observables which are sensitive to the dynamics of 
regions closer to the remnant, where the theoretical situation is not so clear, due to e.g. the 
uncertainties from the BFKL Pomeron and from remnant effects. The observables discussed 
here are defined to probe the current side of the hardest emission in the event. This is a 
relatively clean region with a reduced dependence on soft remnant properties, as compared 
to the forward target region. This facilitates the interpretation of the results, which also can 
be compared to corresponding analyses in e + e~ experiments. We show that the proposed 
observables can distinguish between different existing models. 

The outline of this paper is as follows. In section ^| we discuss suitable jet algorithms in 
a search for high-p^ events. In section the current region properties after different event 
cuts are studied and compared to e + e~ results. In section | we discuss a method to obtain 
experimental e + e~ data whith a controlled rate of heavy flavours. In section | we propose 



observables which are sensitive to the dynamics of regions closer to the remnant. Finally, 
the results of this paper are summarized in section ^| 

2 Jet Algorithms 

In DIS events, high-p^ emissions can break the quark fragmentation universality. In this 
section we discuss jet cluster algorithms suitable for a reconstruction of the highest p_i_-scale 
of an event. The result may depend on the choice of jet finding algorithm used, and we will 
here discuss the properties of a set of possible algorithms. 

Since the aim of the jet clustering is to find the highest p±, it is natural to use k±- 
type cluster algorithms. Our approach will be to combine jets in the order specified by a k± 
distance, until only three remain (including the remnant jet). The event is accepted if the p± 
specified by the final three jets is lower than Q/2. This cut is in analogy to the kinematical 
constraint for gluon emissions in e + e~, which is p± g < E g < a/s/2. 

Historically, k± algorithms were first designed for e + e~ physics. A set of different algo- 
rithms exist In DIS, extra requirements are in general imposed on the cluster algorithm, 
since jet observables in DIS are dependent on structure functions. In a comparison with 
analytical calculations, a jet algorithm where the properties factorizes into perturbatively 
calculable coefficients convoluted with the structure functions is preferred. This implies that 
the remnant jet cannot always be treated on equal footing with the other jets. A k± algo- 
rithm designed to fulfill jet requirements in DIS is presented in Q. However, in the present 
study we are merely interested in jets with relatively high p±, and the infrared properties of 
the used algorithm are less important. Furthermore, we do not intend to investigate jet cross 
sections or the jets themselves, but merely to exclude high-p^ events. In the accepted sample, 
the analysis is performed on the current Breit hemisphere, which is defined independently 
of any jets. We have therefore chosen to use the simpler e + e~ k± algorithms. 

To estimate the reliability of the results obtained using jets, we have used three different 
algorithms, Luclus @, Durham and DlCLUS f|, applied in the hadronic CMS system. 
All these algorithms start off with the final state particles as initial clusters. These clusters 
are then iteratively merged in an order specified by a distance of &j_-type, until all remaining 
distances are larger than a resolution parameter, or until a specified number of clusters 
remain. 

In the HERA experiments, particles which in the lab frame have a large pseudo-rapidity 
w.r.t. the proton direction are not detected. In our analysis, we have chosen to exclude 
all particles with pseudo-rapidity larger than 3.8, to take this into account. For clustering 
purposes, an initial cluster is introduced along the proton direction, carrying the missing 
longitudinal momentum. (To exclude the possibility that the obtained results are merely 
a consequence of the undetected target region, we have also performed analyses where all 
particles are used in the clustering. The result is that the effects presented in the paper 
would be even more prominent, if having access to a full coverage detector.) 

In the Durham algorithm, the distance between the clusters i and j is defined as 



The cluster pair with smallest distance value are combined into one new cluster, with a four- 
momentum given by the sum of its constituents. For small angles 9ij, the Durham distance 





dij is the transverse energy of the softer particle with respect to the other. The factor 
2 sin(0jj/2), used instead of sin(0jj), prevents unreasonable merging of oppositely moving 
clusters. 

The Luclus algorithm is very similar, but uses the distance 

Apart from the difference between |p| and E, the Durham factor min(£ , i , Ej) corresponds 
in the Luclus measure to EiEj/(Ei + Ej). This implies that d^ always is somewhat smaller 

than d\j . To the Luclus algorithm, a reassignment scheme is added, where all particles are 
moved to the closest existing cluster after each merging. At the very end of the algorithm, 
the reassignment procedure is iterated til convergence. This scheme prevents particles to, 
through successive mergings, end up in a jet moving in a very different direction. 

A cluster algorithm can be viewed as a reconstruction of a parton cascade backwards. 
The two described algorithms combines two clusters into one, like an inverse parton emission. 
In the Diclus algorithm, supplied with the Ariadne MC 0, three clusters are merged into 
two. This corresponds to an inverse dipole emission. The distance measure used is 

_ [sij - (rrii + rrij) 2 } [s ik - (m* + m k ) 2 ] 



a i(jk) 



(3) 

Here and are the squared masses of the clusters indicated by the indices. This distance 
is a Lorentz invariant p5_-measure for parton i with respect to a dipole built up by j and k. 
In the merging, the cluster i is absorbed into j and k. The momentum of i is distributed in 
accordance with the recoil treatment in the dipole cascade ||. 

In the present analysis, three clusters are constructed using one of the presented algo- 
rithms. These final clusters are called jets, and the smallest remaining d^ specifies the p± 
of the event. 

The three chosen algorithms are all of k±-tjpe, and hence similar in many respects. 
For quasi-collinear, massless clusters, the three distance measures are very similar. There 
are however also noticeable differences. Differences in the distance measures occur at large 
angles and for heavy clusters, the reassignment procedure in Luclus is not used by the other 
algorithms, and the merging of three clusters into two is a unique feature of the Diclus 
algorithm. In the following, we will investigate results using all three algorithms, and we 
will consider results which are insensitive to the specific choice of algorithm as reliable. 



3 Current Breit Hemisphere Properties 



To examine if a cut in jet-p± is suitable to isolate a sample where universality is expected 
to hold, we will here study results for the current Breit hemisphere in MC generated events 
at HERA energies. The results are compared to e + e~ MC results. 

We simulate the electro-weak interaction in DIS using the Lepto MC ]T0| 
the e + e~ and DIS simulations, we use the Colour Dipole Model 0, 
adne [[|, to describe the parton cascade. The Lund string fragmentation model 
plemented in Jetset |J is then used to describe hadronization. The chosen MC programs 
give a generally good description of data in DIS an e + e~. 



In both 
implemented in Ari- 
im- 
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Figure 1: The rate of empty current Breit hemispheres. Results from MC simulations. The 
solid line shows the rate when considering all generated DIS events. The symbols show the 
result when events with a jet-p± > Q/2 are excluded, using three different k± clustering 
algorithms. The rate is significantly reduced after the cut in jet-p±. 

Effects of the p±-cut 

In Fig |T|, the rate of events in DIS with empty current Breit hemispheres is shown to be 
significantly reduced after a p±< Q/2 cut, for all jet finding algorithms. As discussed in 
section |^, the Luclus distance measure is less restrictive than the Durham one, and this 
is seen as a slightly higher rate of empty events. The rate is however significantly reduced 
with all three algorithms. 

Though the phase space for the current Breit hemisphere is very similar to an e + e~ 
hemisphere, there are differences in the regions nearby, which affect the average energies. 
Consider two hemispheres in the rest frame of an e + e~ annihilation event. In general, they 
have different masses and hence different energies. In other words, a high-p^ emission in one 
hemisphere reduces the energy of the other, due to energy-momentum conservation. 

As for an e + e~ hemisphere, the energy of the current Breit hemisphere at fixed Q 2 is 
expected to deviate from Q/2. The energy depends on emissions in this hemisphere and 
in the nearby phase space of the opposite hemisphere. The latter is in DIS very large, but 
after our suggested event cuts, it is reduced by the condition p±< Q/2. The corresponding 
kinematical constraint in e + e~ is however |p| < Q/2, which is more restrictive. Thus the 
region for high-p^ emissions which can reduce the energy of the considered hemisphere is 
larger in DIS than in e + e~. The emission density at this relatively high p± is suppressed, but 
not zero. Without reaching for a quantitative prediction, we may conclude that we expect 
the mean energy and multiplicity of the current Breit hemisphere to satisfy 

(E CBU ) < |, (iV CB H) < \N ee (Q 2 ), (4) 

where N ee (Q 2 ) is the average multiplicity in an e + e~ experiment with invariant mass Q 2 . 
If 2 (£^cbh) IQ would be much smaller than 1, the similarity between the current Breit 




Figure 2: The relative difference between average energy and Q/2 in the current Breit hemi- 
sphere. The relative shift in the MC simulations is significantly reduced after a p± < Q/2 
cut. 

hemisphere and an unbiased e + e~ hemisphere is poor, and it is not likely that reliable 
conclusions from a comparison can be drawn. On the other hand, if 2 (Ecbu) /Q ~ 1> the 
current Breit hemisphere sample may be closely related to an unbiased e + e~ hemisphere. In 
Fig H, the relative energy shift 2 (-Ecbh) /Q — 1 is shown to be sizeable for the unrestricted 
event sample, but significantly reduced after imposing & p± < Q/2 cut . 

Flavour Compositions 

In a comparison between e + e~ and DIS, also the flavour composition in the different exper- 
iments need to be considered. In e + e~ the rate of heavy flavours varies with energy. On the 
Z° pole it is determined by the couplings to the Z°, while at lower energies it is determined by 
the quark electric charges. In DIS events at HERA energies, photon exchange is dominating 
over Z° exchange. However, here the flavour composition also depends on soft physics. The 
rate of charm in the proton structure function (intrinsic charm) is suppressed by the heavy 
c mass. 

A cut in p± suppresses the boson-gluon fusion channel for charm production in DIS. This 
implies that the c rate in the \ow-p± sample is reduced. Using MC simulations, we have 
found charm to be suppressed approximately by a factor 2 in a p± < Q/2 event sample, 
with all cluster algorithms. A problem for the interpretation of this result is the fact that 
the MC simulation program significantly underestimates the rate of charm events. In the 
MC this is about 10%, while data indicate a much larger value around 25% [13]. However, a 



similar suppression ought to be expected also for real data. If boson-gluon fusion is a more 
important source of charm in Nature than assumed in the MC model, the suppression may 
be even stronger. It would then be appropriate to compare with uds enriched e + e~ data. 
How to obtain such data at different energies will be discussed in section f|. 

Monte Carlo results for average multiplicities are presented in the left plot of Fig [5| 
The results for the low-p^ DIS sample are essentially equivalent for the three considered jet 




Figure 3: MC simulations of multiplicities in e + e~ and the current Breit hemisphere (CBH). 
LEFT: e + e~ results with squared mass s = Q 2 , allowing for all flavours (solid line) and 
only uds (dashed line), compared to CBH multiplicities using all events (crosses) and a jet- 
Pi. < Q/2 sample (diamonds). Both the exclusion of high-p± DIS events and heavy flavour 
e + e~ events improves the agreement between the multiplicities. RIGHT: For low energy 
scales, the agreement is further improved if comparing the CBH multiplicity with e + e~ data 
with squared mass s = 4 (^cbh) 2 - 

algorithms, and only the Diclus result is shown. In this MC generated event sample, heavy 
quark rates are well below 10%, and it is therefore better compared to a MC e + e~ sample 
with only uds. The agreement between the low-p^ and the e + e~ uds sample is better than 
between the unrestricted event samples also shown in Fig |3|. 

Energy Scale Corrections 

As discussed previously in this section, high-p^ emissions close in rapidity to the consid- 
ered hemisphere implies a reduction of the average energy in the current Breit hemisphere, 
(E'cbh), to a value slightly smaller than Q/2. It is then natural to compare the current Breit 
hemisphere multiplicities with e + e~ data at a squared mass s = 4 (.E CB h) 2 - The agreement 
between e + e~ and Breit frame multiplicities is then significantly improved, as shown in the 
right plot of Fig |[ Only the Durham result is presented, but the other results are very 
similar. 

When compensating for small systematic differences by changing scale to the average 
energy, many effects - known as well as unknown - may be corrected for. However, af- 
ter correcting for expected differences in this way, unexpected features which break quark 
fragmentation universality can still be searched for in observables less blunt than the mean 
multiplicity, like strangeness rates, energy spectra and higher multiplicity moments. 

Comparing the left and right plot in Fig |3|, we see that the high-p^ cut has a relatively 
small influence on the average multiplicity, as compared to the energy scale shift. However, 
in Fig [l| we note that the rate of empty current Breit hemisphere events, which are excluded 
from the analysis more by necessity than by theoretical understanding, are reduced with 



the high-px cut. This is also the case for the energy shift, as seen in Fig |2|. Thus we find 
that a cut in jet-pj_ is a powerful step towards an event sample where quark fragmentation 
universality is expected to hold. 

4 Scale Evolutions in Fixed Energy e + e~ Annihilation 

A large sample of uds enriched events are available from LEP1 at y/s = 90GeV. This is not 
the case for other energies corresponding to the HERA kinematical range. In this section we 
discuss how the scale evolution of e + e~ uds hemispheres can be examined, using data from 
a fixed energy e + e~ experiment. 

In search for an algorithm which enables such a study, it is useful to recollect the Breit 
frame analysis discussed previously in this paper. In that analysis three jets are found with 
a k±_ algorithm in the hadronic center of mass frame. In events were the p± of the event is 
smaller than Q/2, the current Breit hemisphere is investigated. This region corresponds to a 
cone in the hadronic center of mass frame, defined by a constraint on the rapidity measured 
in the current direction, 

y>~ln(W 2 /Q 2 ), (5) 

where W 2 is the squared mass of the hadronic system. 

We will now consider a very similar algorithm for e + e~ . In an e + e~ experiment at 
squared mass s, we consider an artificial scale Q max < s. Three jets are reconstructed with 
a k± cluster algorithm, and events where p±> Q max /2 are excluded. In the accepted event 
sample, Q maiX /2 can be identified with the maximal allowed virtuality for a parton in the 
event. Thus the multiplicity depend on two scales, energy and virtuality. For Q max = s, the 
two scales coincide and we have an unbiased event sample. 

We study the multiplicity for particles where the rapidity, measured in the thrust direc- 
tion, satisfies 

\y\ > hn(s/Q 2 m J. (6) 

In e + e~, where no proton remnant is present, both the positive and negative rapidity cone 
is considered. These cones in the thrust direction correspond to e + e~ hemispheres with 
squared invariant mass Q max , and their evolution with this scale can thus be studied in a 
fixed energy experiment. (For a more detailed discussion, see e.g. 

The "thrust-cone" method presented here is designed to be as similar as possible to the 
Breit frame analysis. In both e + e~ and DIS, jets reconstructed in the hadronic rest frame are 
used to exclude high-p^ events. In the accepted events, the analysis is performed in rapidity 
cones defined by Eq fl5|) and Eq (||), respectively. These cones correspond in some other 
Lorentz frame to unbiased hemispheres, where the maximal allowed transverse momentum 
coincides with the hemisphere energy. (In DIS, the boost to the relevant Lorentz frame, the 
Breit frame, is performed explicitly.) 

Two differences are present due to the absence of a t-channel probe in e + e~ . The scale 
Q 2 and the direction which defines the rapidity cones are determined by the probe in DIS. In 
e + e~, the scale Q max is chosen freely and the rapidity is measured along the thrust direction. 

The upper limit on p]_ in a jet corresponds to the maximal allowed virtuality for the 
parton initiating the jet. Jet definitions suitable for a systematic study of the dependence 




Figure 4: Multiplicities in MC simulated e + e~ events including all allowed flavours (dotted 
line) are higher than for light quarks events (dashed line). The multiplicity obtained by the 
thrust-cone algorithm in a uds sample at the Z° pole (solid line) are in better agreement 
with the uds result for moderate scales between 4 and 8 GeV. This corresponds to current 
Breit hemispheres in DIS with Q 2 between 15 and 60 GeV . For larger scales the thrust-cone 
results differ from the uds sample in a similar way as the full flavour results. 



upon the two separate scales, energy and virtuality, are presented in |TJ]], and will be further 
discussed in the section || of this paper. In the comparison of e + e~ results with data from 
the current Breit hemisphere in DIS, we are however mainly interested in quark jets where 
the two scales coincide. It is then sufficient to study the scale evolution using the relatively 
simple thrust-cone method. 

In Fig |] the multiplicity evolution of an e + e~ uds sample is compared to results including 
all flavours and results obtained by the thrust-cone algorithm in a fixed energy uds sample. 
The thrust-cone multiplicities are plotted as a function of the average energy scales obtained 
in the cones. These scales are somewhat smaller than the chosen Q max /2, and the deviations 
are of order 5 to 10%. 

In rare cases where an emitted gluon gives the hardest jet, the thrust will be directed 
along the gluon momentum. Thus the average multiplicity in the thrust-cone analysis will 
get some contribution from gluon jets, which have a larger multiplicity. This results in a 
systematic overestimation of the multiplicity, as shown in Fig [|. 

The aim of the thrust-cone algorithm is to reproduce the scale evolution of a e + e~ uds 
sample better than a sample including all flavours. As seen in Fig [|, this is achieved for 
moderate Q max , between 4 and 8 GeV. This corresponds to current Breit hemispheres in DIS 
with Q 2 between 15 and 60 GeV 2 , which is an important range at the HERA experiments. 

For Q max larger than 8 GeV, where large p± gluon jets are not excluded, the effect 
of gluon "pollution" is of the same order as the effect of the heavy quarks. However, we 
note that Fig |] shows results from MC simulations. In a realistic experimental analysis, 
the thrust-cone results applied on uds events at the Z° peak will benefit from very large 
statistics compared to experiments at lower energies, and also from the fact that the scale 



evolution can be studied over a large range with the same detector. To conclude this section, 
our investigation indicates that it would be interesting to compare DIS data not only to full 
e + e~ results at different s, but also with thrust-cone results in uds samples from the Z° 
experiments. 



5 Approaching the Proton Remnant 

Compared to the current Breit hemisphere, the target region in the Breit frame is theo- 
retically much less understood. It is difficult to determine the basic dynamics (e.g. the 
importance of the BFKL mechanism) using only inclusive observables, like i^-flow. In the 
far forward direction of the target region, results depend to a large extent on non-trivial soft 
properties of the remnant. This region is also more difficult to study in many experiments. 
E.g. at HERA, many particles are lost in the beam pipe. 

In this section we discuss two observables which examine the current and central rapidity 
ranges, but which nevertheless are sensitive to soft remnant effects and the dynamics of the 
evolution. The first observable is actually the same one that has been discussed throughout 
the paper, i.e. the multiplicity in the current Breit hemisphere. We will however here focus 
on results for very high Q 2 . The rapidity distance to the remnant is then reduced, and as 
we will see, different models for remnant effects predict different results in this observable. 

We will also propose a method to gradually extend the investigated phase space region, 
from the current Breit hemisphere, to also include central rapidity regions of the target 
hemisphere. Compared to the forward target region, this region is less dependent on soft 
remnant effects, but the basic mechanism for radiation is not very well known and depends to 
some extent on rather subtle coherence effects (cf. the "Feynman-Gribov puzzle" discussed 



in pi). 



We will use a jet definition, "the Mercedes algorithm" |1J], with which we can construct 
quark jets with well-defined rapidity regions both in DIS and e + e~ . This algorithm can be 
used to examine current jets in DIS, where the considered rapidity range is large enough 
to reach into the target hemisphere of the Breit frame, and the result can be compared to 
similarly defined jets in e + e~ . 

The Mercedes algorithm in DIS is applied to three-jet configurations, with a current jet, a 
remnant jet and a third jet. Around the current jet a cone-like region is studied, which covers 
rapidities up to the third resolved jet. If this is found outside the region which defines the 
current Breit hemisphere, as illustrated in Fig [5], the Mercedes current jet will probe central 
rapidities in the target hemisphere of the Breit frame. It will however only cover rapidities 
on the current side of the highest p± emission. This is a comparatively clean region in which 
relatively few dynamical mechanisms are important, as compared to more forward regions 
of the target hemisphere. 

We will in this section start to present the Mercedes algorithm, applied e + e~ and DIS, 
respectively. The scale dependences and expected properties of the Mercedes jets are then 
discussed. After a presentation of some different models in DIS, we finally show that Mer- 
cedes jet results, and also the current Breit hemisphere results at high Q 2 , can be used to 
distinguish between these. 
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Figure 5: The current Breit hemisphere roughly corresponds to a cone around the current 
jet in the hadronic center of mass frame (light-gray cone). With the Mercedes algorithm, a 
current jet region defined to cover the rapidity range up to a resolved high-p± jet is studied 
(both gray cones). In the figure the third jet lies outside the current Breit hemisphere region 
and the Mercedes current jet will in this example include central rapidities in the target 
hemisphere of the Breit frame. 

5.1 The Mercedes Algorithm 

The starting point of the Mercedes algorithm is a three-jet configuration in e + e~, obtained 
by a k± cluster algorithm in the center of mass frame. The p± of the e + e~ event is specified as 
the smallest cluster /c_i_- distance between the remaining jets. This implies that no emissions 
within the jets have higher transverse momentum, and p\ can thus be identified with the 
maximal virtuality of a parton initiating a jet in the event. 

The jets are analysed after a boost to the symmetric "Mercedes" frame, where all jet 
angles are 9' = 120°. As the jet definition is not Lorentz invariant, such a Lorentz trans- 
formation can transfer soft particles from one jet to another. Thus the mass of a jet is not 
invariant. Instead the jet direction can be interpreted as the direction of a massless parton 
initiating the jet, which implies that it transforms as a light-like vector. (The Diclus algo- 
rithm is a Lorentz-invariant cluster algorithm, which thus is an exception to this discussion. 
However, particles are not assigned to any specific jet in this algorithm, and the produced 
jets are explicitly massless.) 

In the Mercedes frame the boundaries between the jets are determined by the bisectors 
between them. In [JTJJ it is shown that when the gluon is significantly softer than the quarks, 
the rapidity ranges (at a lower cut-off scale A) for the quark jets are given by 

(7) 

The first term on the right-hand side of Eq (|?p is the rapidity range in the forward region 
of the jet where 2y > ln(s/4pj_), which corresponds to an unbiased hemisphere of an event 
with squared invariant mass equal to Ap\ (cf. Eq (El)). The second term is the extra rapidity 




range included in the Mercedes quark jet. In this range the maximal allowed transverse 
momentum is p±. 

Thus the jet properties depend on two independent scales. The virtuality scale p± de- 
termines the size of the first range and the maximal allowed transverse momentum in the 
second, while the energy sets the size of the second range. This is not a unique feature of the 
Mercedes jet, all jets can in general be subdivided into two regions in a similar way, and thus 
depend on two scales. When looking at a forward cone, so that the second term vanishes, 
the two scales coincide. The cone then corresponds to an unbiased e + e~ hemisphere. In the 
region corresponding to the second term in Eq (0), where the maximal transverse momentum 
is constant, the emission density in an e + e~ jet is essentially flat in rapidity. This implies 
that the multiplicity in the jet is expected to grow linearly with the allowed rapidity range, 
for fixed virtualities. 

When applying the Mercedes algorithm to DIS events, we will use the Durham k± cluster 
algorithm in the hadronic center of mass frame, in order to resemble the corresponding e + e~ 
analysis. Three jets are reconstructed, and the smallest remaining /cj_-distance between these 
specifies the p± of the event. The current jet is tagged as the one with smallest polar angle 
to the probe direction, and is analyzed in the Mercedes frame. The considered rapidity range 
for the current jet is then given by Eq (|7]), with x q being the scaled energy 2E q /W of the 
outgoing struck quark in the hadronic center of mass frame. 

The Mercedes current jet can reach into the target region of the Breit frame. The 
rapidity range for the current Breit hemisphere (the light-gray cone in Fig ||) is \n(Q/A). 
Eq (|7D implies that the Mercedes current jet in DIS will cover a region which extends into 
the target hemisphere of the Breit frame (the dark-gray cone in Fig [5]) when p± > Q(l —x q ). 

As for the e + e~ case, the properties of the Mercedes current jet in DIS depends both 
on its energy and its virtuality. It may, however, also depend on the DIS scales x and Q 2 . 
These scale dependencies can be investigated for a set of observables. E.g., the Mercedes 
jet properties in events where p±> Q could be used to investigate the mechanisms behind 
resolved photon events (16|. It is also possible to investigate different rapidity intervals in 
the Mercedes current jet separately. 

Here we will investigate whether the Mercedes jet can be used to distinguish between 
different models in DIS. It is then sufficient to study the simplest possible observable, the 
average multiplicity in the whole jet. We will compare the rapidity range dependence of 
the average multiplicity in the DIS current jet and e + e~ quark jets, keeping the virtuality 
scale fixed. We note that when Q is smaller than p±, a significant Q 2 dependence could be 
expected for the current jet in DIS. For the jets in the two experiments to be similar, we 
therefore exclude DIS events where p± > Q/2. In the remaining event sample, we study the 
inclusive results, integrated over Q 2 and x. 



5.2 Models 



To demonstrate that the proposed observable can distinguish between different production 
mechanisms, we study some different models and see that they indeed give different results. 
We have chosen two versions of the Lepto MC |TD|, based on a DGLAP parton evolution, 
and two versions of the Ariadne MC 0, implementing the soft radiation model [17j with 
boson-gluon fusion included. All the chosen models uses the Lund string fragmentation to 



describe hadronization. The models however differ in the treatment of parton evolution and 
proton remnant effects. 

In the Lepto MC, a parton shower based on the DGLAP evolution equation is imple- 
mented. This implies that parton emissions are strongly ordered in rapidity and transverse 
momentum. A ladder diagram describing a DIS event specifies a colour ordering for the 
emissions, which due to colour coherence also specifies the order in rapidity. For large Q 2 , 
the leading order contribution to the cross section is given by diagrams where the transverse 
momenta are strongly ordered, so that the parton closest in rapidity to the remnant has the 
smallest p±, and the parton closest to the quark interacting with the probe the largest. It 
is diagrams with this ordering in transverse momentum that are considered in the DGLAP 
evolution. 

A model for soft colour interaction (SCI) is supplied in the Lepto MC. In this model, 
the proton remnant may interact with the perturbatively produced partons by exchanging 
soft colour quanta. These have negligible effect on the parton momenta, but alter the colour 
flow of the parton configuration. This implies that a qq pair produced in a boson-gluon fusion 
event may change its colour state from an octet to a singlet. This introduces a rapidity gap in 
the colour flow from the current to the remnant. In cluster and string fragmentation models 
the colour topology determines the phase space for hadron production, and thus the same 
rapidity gap manifests itself in the hadronic final state. This is a way to model diffractive 
DIS events, alternative to the Pomeron approach. 

Another possible consequence of SCI is that the colour flow is enlarged, thus producing 
more hadrons. This increases the E^-flow in a way which is consistent with data |18| . 

In the default version of the Ariadne MC, the radiation in DIS is expressed in terms 
of a dipole spanned from the struck quark to the remnant. The emitted partons are not 
necessarily ordered in k± and rapidity as in the DGLAP case, but the fact that the remnant 
is an extended object, and not point-like, gives rise to a suppression of high-p^ emissions in 
the remnant direction. Though this model is rather different from the DGLAP formalism, 
the radiation from the struck quark in the current region is treated similarly in the two 
Monte Carlos. 

The suppression of high-p^ emissions, in the way it is modeled in Ariadne, may reach 
into the current Breit hemisphere when Q 2 is very high. It could however be argued that 
the suppression in this case is overestimated, since the maximal allowed virtuality for the 
struck quark then is lower than Q 2 . A modification has therefore been developed, where 
the unsuppressed phase space is enlarged in high Q 2 events, to always include the expected 
region for final state emissions from the struck quark [IK . 



5.3 Results 

We will here study the prediction of different MC simulation models on two observables, 
the multiplicity in the current Breit hemisphere at high Q 2 and the current jet multiplicity 
obtained by the Mercedes algorithm. 

In Fig the prediction for the multiplicities in the current Breit hemisphere are shown, 
using different models. At very high Q 2 , the rapidity distance to the remnant is reduced, 
and different assumptions about remnant effects give different predictions . 

The differences can be understood from the remnant treatments of the models. The 
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Figure 6: Different model predictions for the multiplicity in the current Breit hemisphere. In 
the Ariadne default version the remnant suppression of high-p^ emissions can reach into 
the current Breit hemisphere in events with high Q 2 . This implies a lower multiplicity. In 
the Lepto MC with soft colour interaction, where the hadronization phase space in terms 
of a colour string can be enlarged, the multiplicities lie above the others. 

Ariadne default version, where the remnant suppression of high-p^ emissions can reach 
into the current Breit hemisphere for very large Q, predicts a lower multiplicity than the 
other models. The prediction from Lepto with soft colour interactions (SCI) lies above the 
others. SCI can give rise both to rapidity gaps and an extended colour flow, which increases 
the hadron multiplicity. The results here indicates that the latter effect dominates in the 
current rapidity region. 

The Lepto MC without SCI and the high-Q modified Ariadne MC differ mainly in 
the evolution of the parton distributions. In Lepto the evolution is based on the DGLAP 
equation, while the Ariadne MC implements a soft radiation model based on colour dipoles. 
The properties of the current Breit hemisphere are expected to be essentially independent of 
the form of the parton distribution evolution. The two MC results are similar at very high 
Q 2 , in agreement with this expectation. 

The MC result for an e + e~ uds sample (not shown) agrees well with the high-Q modified 
Ariadne and Lepto without SCI. Whether this agreement between e + e~ and DIS results 
with similar flavour compositions holds also for experimental data is however not fully inves- 
tigated yet. The differences between the results in Fig |6] show that the high-Q 2 Breit frame 
analysis can be used to discriminate between different models. 

Fig shows the Mercedes algorithm multiplicity of the current jet in DIS and quark 
jets in e + e~, as a function of the available rapidity range y. The transverse momentum is 
constrained to the range 3GeV< p± <3.6GeV and DIS events with Q/2 < p±_ are excluded 
from the analysis. We have investigated the results for a rather low p±, in order to make 
the analysis complementary to the high-Q 2 analysis shown in Fig The multiplicity of the 
e + e _ jets rise fairly linearly with y, as expected. For low values of y the examined region 
lies within the current Breit hemisphere, and the DIS results are similar to the e + e~ result. 
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Figure 7: MC results on the multiplicity in DIS current jets and e + e~ quark jets using the 
Mercedes algorithm, as a function of the rapidity range in Eq (ffi). The transverse momentum 
scale is constrained to 3 GeV < p±_ < 3.6 GeV. 

Ideally, all jets correspond to a perturbative gluon emission, with similar p±. However, 
"junk jets" may be formed, where hadrons which happen to be close in phase space are 
merged into a resolved jet, even though no corresponding gluon was emitted. These "junk 
jets" are more easily formed when the multiplicity is high. Thus their presence in the 
considered event sample increases the average multiplicity. If the probability to emit gluons 
giving rise to "proper" jets is reduced, the fraction of junk jet events in the accepted sample 
increases, and their effect on the result gets more prominent. The results for Mercedes jets 
with large rapidity ranges are based on events where one jet is found in the target region. 
There the probability for gluon emission is smaller than in a corresponding region in e + e~, 
and in Fig [7] the junk jet contribution is seen as a larger multiplicity for DIS Mercedes jets 
with large rapidity ranges. 

The probability for jet-generating gluons in the target region differ in the considered DIS 
models, being lower for the DGLAP based Lepto MC than the dipole based Ariadne MC. 
Thus the high-Q modified Ariadne MC and the Lepto MC without SCI, which give very 
similar results in the high-Q 2 Breit frame analysis of Fig |^, differ in the Mercedes algorithm 
results in Fig ^. This result suggests that properties of the parton evolution, such as p± 
ordering, can be examined using the Mercedes algorithm for jet analyses. 



6 Summary 

The assumption of quark fragmentation universality implies that the current Breit hemi- 
sphere in DIS is expected to be very similar to a hemisphere in e + e~ annihilation. However, 
the experimental situations are different, and several corrections to universality are present. 
We have here proposed event cuts to improve the expected validity of quark fragmentation 
universality. DIS events with high-p^ emissions, which have no correspondence in e + e~ 
events, are excluded, where the p_i_-scale is reconstructed using jet cluster algorithms. 



The effects of the cuts have been investigated using MC simulations, and the agreement 
between e + e~ and Breit frame results is shown to improve. The results after a cut in jet p± 
may depend on the choice of cluster algorithm. We have investigated three different types 
of k± cluster schemes, and find our results to be algorithm independent. 

In the accepted low-p^ DIS sample, heavy quarks are suppressed. This motivates a 
comparison with uds enriched e + e~ data, which are available from the experiments at the 
Z° pole, but not at lower energies. We have here presented a method, the "thrust-cone 
algorithm", to study scale evolutions of e + e~ quark hemispheres, using data from fixed 
energy experiments. With this algorithm, uds enriched data with high statistics from the 
LEP1 experiments can be compared to results for the current Breit hemisphere, over a large 
range of energies. 

At high Q 2 , the rapidity distance from the current Breit hemisphere to the remnant end 
is reduced. By comparing results from different MC programs, we find that the multiplicity 
in the current Breit hemisphere at very high Q 2 can be used to distinguish between different 
assumptions of remnant effects. 

We also propose an observable which more directly investigates the properties of the 
target region. A current jet is defined with the "Mercedes" algorithm [|14]]. The Mercedes 
current jet probes central rapidity regions in the Breit frame, but always stays on the cur- 
rent side of the highest p± emission of the event. This is a comparatively clean region, as 
compared to the forward target region. A corresponding jet definition can be studied in 
e + e~ experiments, which enables a comparison between the processes, and could facilitate 
the interpretation of the results. The properties of the Mercedes current jet is shown to be 
sensitive both to remnant effects and mechanisms behind the parton evolution, and can be 
used to discriminate between different existing models. 
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